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Resonant interaction of a collection of two-level atoms with 
a single-mode coherent cavity field is considered in the frame- 
work of the Dicke model. We focus on the role of collective 
atomic effects in the phenomenon of collapses and revivals of 
the Rabi oscillations. It is shown that the behavior of the 
system strongly depends on the initial atomic state. In the 
case of the initial half-excited Dicke state we account for a 
number of interesting phenomena. The correlations between 
the atoms result in a suppression of the revival amplitude, and 
the revival time is halved, compared to the uncorrelated fully- 
excited and ground states. The phenomenon of squeezing of 
the radiation field in the atom-field interaction is also dis- 
cussed. For the initial fully-excited and ground atomic states, 
the field is squeezed on the short-time scale, and squeezing can 
be enhanced by increasing the number of atoms. Some em- 
pirical formulas are found which describe the behavior of the 
system in excellent agreement with numerical results. For the 
half-excited Dicke state, the field can be strongly squeezed on 
the long-time scale in the case of two atoms. This kind of 
squeezing is enhanced by increasing the intensity of the ini- 
tial coherent field and is of the same nature as revival-time 
squeezing in the Jaynes-Cummings model. The appearance 
of this long-time squeezing can be explained using the factor- 
ization approximation for semiclassical atomic states. 
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I. INTRODUCTION 

Since the pioneering work of Dicke [Q on cooperative 
spontaneous emission, a great deal of attention has been 
devoted to the interaction of the radiation field with a 
collection of two-level atoms located within a distance 
much smaller than the wavelength of the radiation. Such 
a system is commonly referred to as the Dicke model 
(for a review see, e.g., Ref. A particular case of the 
Dicke model, when atoms interact with a single-mode 
radiation field inside a cavity, was considered by Tavis 
and Cummings J3|. The Tavis- Cummings Hamiltonian is 
mathematically equivalent to the trilinear boson Hamil- 
tonian describing various nonlinear optical processes ||. 
The single-atom version, known as the Jaynes-Cummings 
model ||, is the simplest and one of the most popular 
models of quantum optics. In spite of its simplicity, the 
Jaynes-Cummings model shows a variety of interesting 
nonclassical phenomena such as vacuum-field Rabi oscil- 
lations, sub-Poissonian photon statistics, and squeezing 
of the radiation field (for reviews see, e.g., Refs. 

One of the most interesting quantum features of the 
Jaynes-Cummings model is the phenomenon of collapses 



and revivals of the Rabi oscillations, which manifests it- 
self in the clearest way when the cavity field is prepared 
initially in the coherent state ||. The shape of collapses 
and revivals is determined by the initial photon-number 
distribution. A similar behavior can be found also in the 
many-atom case [p|-p|. For a sufficiently strong coher- 
ent field, the nonlincarity in the Rabi frequency is slight 
and the system exhibits regular dynamics in the form of 
collapses and revivals of the oscillations. However, in the 
many-atom case there exist anharmonic collective cor- 
rections which modify the shape of the collapses and re- 
vivals related to the photon-distribution mechanism [fl2| . 
In the present work we study in detail these collective 
effects for different initial atomic states. If atoms are 
prepared initially in the ground state or in the fully- 
excited state, then the system behaves rather similarly 
to the single-atom case, although collective effects man- 
ifest themselves clearly in the length of the revival time. 
But, the results are different when the atoms are pre- 
pared initially in the half-excited Dicke state. This state 
is well known as the superradiant atomic state in the 
context of collective spontaneous emission in open space 
It was also found (r| that the half-excited Dicke 
state strongly modifies the resonance fluorescence spec- 
trum of the atomic system. However, to the best of our 
knowledge, this state was rarely considered in the con- 
text of the Tavis-Cummings (cavity) version of the Dicke 
model. When the atoms are prepared in the half-excited 
Dicke state the system exhibits a number of interesting 
phenomena. In particular, the amplitude of the Rabi os- 
cillations is strongly suppressed and relative intensities of 
revivals are essentially changed. This suppresion of the 
revival amplitude is similar to the trapping phenomenon 
occuring in the case of a single atom prepared in the 
equally weighted superposition of the two levels jyj. 

It has been known for a long time fn|[L(| that the non- 
linear character of the Jaynes-Cummings model leads to 
squeezing in one of the quadratures of an initially coher- 
ent cavity field. It was also predicted jlTj that strong 
squeezing can be obtained in the Jaynes-Cummings 
model near the revival times for large initial intensities 
of the field. This phenomenon of strong revival-time 
squeezing was explained in Ref. [ fl8[ using the factoriza- 
tion approximation for the so-called semiclassical atomic 
states fl9f| . Butler and Drummond |^(| showed that 
short-time squeezing can be enhanced in the Dicke model 
compared to the single-atom case. Seke considered 
field and atomic squeezing in the Dicke model without 
the rotating-wave approximation and in the presence of 
losses. Analytical approaches to squeezing in the strong- 
field limit and in the weak-field limit ||^] were pre- 
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sented. Higher-order squeezing in the Dicke model was 
studied in Ref. H]. For atoms prepared initially in the 
ground state or in the fully-excited state, the uncertainty 
of the field quadrature rapidly oscillates and squeezing 
occurs during short periods of time at the very beginning 
of the time evolution. On the long-time scale the tempo- 
ral behavior of the quadrature uncertainty is correlated 
with times of collapses and revivals of the Rabi oscil- 
lations. On the other hand, for the initial half-excited 
Dicke state, the behavior of squeezing is qualitatively 
different. The quadrature uncertainty oscillates on the 
long-time scale, with a period of the order of the revival 
half-time, and strong squeezing is obtained in the case of 
two atoms. The value of squeezing is enhanced by the 
intensity of the initial coherent field. This phenomenon 
is of the same nature as revival-time squeezing in the 
Jaynes-Cummings model and can be explained using the 
factorization approximation for the semiclassical atomic 
states. In general, squeezing is better and its duration is 
much longer for the initial two-atom half-excited Dicke 
state than for the ground and fully-excited states. 



II. THE MODEL AND METHODS OF SOLUTION 

We consider the resonant interaction between N two- 
level atoms and the single-mode radiation field inside a 
lossless cavity. In the rotating-wave approximation, the 
Tavis-Cummings interaction Hamiltonian reads {fi = 1): 



H = gia^J- + aJ + ). 



(2.1) 



Here an d in the following we use the interaction picture. 
In Eq. (|]l]) g is the coupling constant, a and are the 
annihilation and creation operators of the field mode, J+ 
and J- are the collective atomic raising and lowering 
operators. They satisfy the su(2) Lie algebra, 



[J+, J-]=2J Z 



[J„J±]=±J±, 



(2.2) 



where J z is the operator of atomic inversion. In terms 
of the standard Pauli matrices, describing each two-level 
atom, one obtains 



J± 



1 N 



1 



A" 



i=l 



(2.3) 



If the Schwinger realiza tion of the SU(2) generators is 
used, the Hamiltonian (2.1) becomes the trilinear bo- 
son Hamiltonian describing nonlinear optical processes 
such as parametric conversion and Raman and Brillouin 
scattering Q . An interesting physical realization of the 
Hamiltonian (2.1) is given by the coupling of the inter- 
nal levels of atoms or ions to a mode of their quantized 
oscillatory motion in a harmonic trap p5[ . 
The total excitation operator 



commutes with the Hamiltonian (2.1) and is an inte- 
gral of motion. Another integral of motion is the SU(2) 
Casimir operator J 2 = J 2 + \{J + J- + J_ J+). We de- 
scribe the state of the atomic system in terms of the 
SU(2) orthonormal basis \j, m) a t (to = j,j — 1, . . . , — j), 



Jz\j,m) at = TO | J, TO,) a t, 
J 2 |j,TO) at =j(i + l)|j,TO) at . 



(2-5) 
(2.6) 



In the context of the atomic realization (2.3), the states 
|j, to) a t are the symmetric Dicke states: 



1 3, rn) 



-1/2 



e n i+>«. n 



(2.7) 



k=l 



where |+); and |— ); are the upper and lower states, re- 
spectively, of the Zth atom, and the summation is over 
all possible permutations of N atoms. If only symmetric 
atomic states are considered, then the 'cooperative num- 
ber' j is equal to N/2 and p = m + j is just the number 
of excited atoms. 

The Hilbert space TC of the atom-field system can be 
decomposed into a direct sum of finite-dimensional in- 
variant subspaces TIl- 



L = a)a + J z + N/2 



(2.4) 



L=0 

Each invariant subspace TIl is spanned by the orthonor- 
mal basis \n)f\j,L — j — n) a ti where \n)f are the Fock 
states of the radiation field, a<a\n){ = n\n)f. For L < N, 
n = 0,1,..., L and dim(H L ) = L + 1; for L > N, 
n = L- N, L-N+1,...,L and &mx(TL L ) = jV+1. If the 
field is initially in the Fock state |no)f and the atoms are 
in the Dicke state \j, mo) a t, the state of the system will 
evolve in the invariant subspace TIl with L = no+TOo+j. 
For the field and/or atoms prepared initially in a super- 
position state, one should take into account contributions 
from different subspaces. 

The exact solution of the problem is obtain ed b y the di- 
agonalization of the interaction Hamiltonian ( |2.l| ) in each 
of the invariant subspaces TLl involved |3| . It is known 
|| that in the basis \n) { \j,L — j — n)a,t the Hamiltonian 
is given by a tridiagonal matrix with symmetric eigenval- 
ues and the corresponding characteristic equation can be 
reduced to an algebraic equation of order [dim(W/,)/2]. 
Therefore, an analytical solution is possible when only in- 
variant subspaces with dim(7ii) < 9 are involved. How- 
ever, already for dim (7^) > 3, analytical solutions are 
rather complicated [Pq-p8|. Note that for dim(H.L) > 4 
the eigenvalues are not equidistant, so the time evolution 
is not periodic even within a single subspace. Semiclas- 
sical approximate solutions were proposed p9Lpl[ which 
give the time evolution of systems governed by a trilinear 
Hamiltonian of type (|2.1| ) in terms of elliptical functions. 
A perturbative analytical approach to the problem with 
weak fields was developed by Kozierowski and co-workers 
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p7| , p8 32 3j|. In the present work we use the exact solu- 
tion based on the nu mer ical diagonalization of the inter- 
action Hamiltonian (2.1). 



III. COLLAPSES AND REVIVALS 

We study the temporal behavior of the atom-field sys- 
tem in the Dicke model for the cavity field prepared ini- 
tially in the coherent state \a)f. 



ah 
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71=0 



n){. 



(3.1) 



Without loss of generality we consider only real values of 
a. The initial atomic state is supposed to be one of the 
Dicke states \j, m) a t (recall that j = N/2). Two possi- 
bilities which are frequently considered in the literature 
are the fully-excited state |j, j) a t and the ground state 
\j, —j) a t- We are also interested in the half-excited Dicke 
state \ j, 0)at- 

Since the coherent field state is a superposition of many 
Fock states \n){, the invariant subspaces Hl with differ- 
ent values of L = n + m + j contribute to the evolution. 
The temporal behavior of the atomic inversion (J z ) is 
given by the sum of the appropriately weighted atomic 
responses to each Fock state. (The mean photon num- 
ber (a^a) is connected with the atomic inversion by the 
fact that (L) = (a T a) + (J z ) + N/2 is a constant.) In 
the resonant Jaynes-Cummings model (N — 1) with ini- 
tially unexcited (to = — i) or excited (to — -|) atom, one 
obtains lp| 



{Jz 



■ ^ y Pfi COS £i n rn Tj 



where t = gt is the scaled time, 

&n,m = 2y/n + TO 



1/2 



(3.2) 



(3.3) 



is the scaled Rabi frequency corresponding to each sub- 
space, and P n is the photon-number distribution. For the 
initial coherent state, P n = exp(— n)n n /n\ is the Poisso- 



nian distribution and n 



is the initial mean pho- 



ton number. Then, due to the property of this distribu- 
tion, the most regular dynamics occurs for large values 
of the initial mean photon number n. Contributions cor- 
responding to different n's interfere in such a manner 
that they initially go out of phase, after that acquire a 
common phase, and this process is repeated resulting in 
a series of collapses and revivals, as shown in Fig. [y(a). 
The revival time tr can be estimated using the condition 



n+l,m 



^n,m) — ^7T, 



(3.4) 



which gives tr = 2it\J n + m + 1 /2. 



One should expect a similar behavior also for the Dicke 
model in the strong- field domain n > N - 12 1 . Of 
course, the collapses and revivals related to the photon- 
number distribution would be modified by the collec- 
tive atomic effects due to the fact that the eigenvalues 
of the interaction Hamiltonian are not equidistant. In 
the strong-field limit n 3> N the anharmonic correc- 
tions to the eigenvalues become small, and one should 
expect a quite regular behavior similar to that of the 
resonant Jaynes-Cummings model. However, we will see 
that there exist reasons for an irregular behavior that 
are much more important than just the anharmonicity of 
the eigenvalues. In fact, the initial atomic state deter- 
mines how important will be various factors leading to 
irregularities in the behavior of the system. 
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FIG. 1. The atomic inversion (J z ) versus the scaled time 
t — gt for the initial coherent field state with n = 36 in- 
teracting with (a) one excited atom (the Jaynes-Cummings 
model) and with 10 atoms prepared in (b) the fully-excited 
state \j,j)ax and (c) the ground state \ j, — j) a t. 

We first consider the cases with the atoms prepared 
initially in the fully-excited and ground states. The tem- 
poral behavior of the atomic inversion (J z ) is shown in 
Fig. [l](b,c) for the case of N = 10 and n = 36, with atoms 
prepared in (b) the fully-excited state \j,j) a t an d (c) the 
ground state \j, — j) a t- The first feature to note is the de- 
pendence of the revival time on the initial atomic state. 
Similarly to the Jaynes-Cummings model, we can esti- 
mate the revival time for the initial atomic state \ j, m) a t 
as 



2-K^ff, 



1/2, 



(3.5) 



where we use the strong-field limit expression for the Rabi 
frequency (i.e., neglect the anharmonic corrections to the 
eigenvalues). Formula ( |3.5[ ) is in a good agreement with 
our numerical results, as demonstrated in Fig. ^[ The 
difference between the revival times for the fully-excited 
state and the ground state is particularly obvious when 
N is not too small compared to n. The second feature is 
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that the spread of the revivals is determined only by the 
photon statistics of the initial field state 34 and does not 
depend on the value of to. These two observations may 
explain why the temporal behavior loses its regularity 
for the initial ground state notably faster than for the 
initial fully-excited state. Indeed, we see in Fig. [j](b) 
four regular revivals for m — j, while in Fig. UKc) for 
to = —j the behavior starts to be irregular already at the 
end of the third revival. The reason is that the spread of 
the revivals increases with time and neighboring revivals 
start to overlap, which leads to the loss of regularity. 
Since for m = —j the revival time is shorter than for to = 
j, the overlapping of the revivals and the corresponding 
irregularity occur earlier. 
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FIG. 2. The atomic inversion (J z ) versus the scaled time 
t — gt for the initial coherent field state with n = 36 inter- 
acting with with 10 atoms prepared in (a) the fully-excited 
state \j,j)at, (b) the half-excited Dicke state \j, 0) a t, and (c) 
the ground state \j, — j)at. The region of the first revival is 
shown, demonstrating the dependence of the revival time on 
the initial atomic state. 

Another factor that leads to differences between the 
behavior of the system in the cases to = j and to = — j 
is the dependence of the anharmonic corrections to the 
eigenvalues on the value of to. In order to give a rep- 
resentative example, we considered the eigenvalues hi 
of the interaction Hamiltonian for the subspaces with 
L = n + m + j for n = 36, N = 10, 20, 30 and to = ±j, 0. 
The ratio /51 = hs/hi of the 5th and 1st eigenvalues 
was chosen as a parameter representing the anharmonic- 
ity (for equidistant eigenvalues /51 is exactly 5). The 
numerical results are listed in Table |. We see that for 
given N the value of /51 increases with decrease of to, i.e., 
irregular effects related to the anharmonicity are most 
important for to = —j. For to = j the value of /51 de- 
creases slightly with increase of N, but the overall con- 
tribution of the anharmonic corrections increases with N 
just because there are more eigenvalues. For to = — j 
the value of /51 increases with AT", so the importance of 
the anharmonic corrections here increases with N much 



faster than for m = j. In particular, we see that the 
amplitude of the revivals for to = — j is smaller than for 
to = j, and this effect becomes more pronounced as N 
increases. For n = 36, as N increases from 2 to 16, the 
relative amplitude Ai (the difference between the max- 
imum and minimum values of (J 2 ) in the first revival, 
divided by N) decreases from 0.5459 to 0.4039 for m — j 
and from 0.5362 to 0.2112 for to = — j. This effect occurs 
because the anharmonic corrections partially destroy the 
interference of the oscillating terms. Of course, as the 
initial mean photon number n increases, the behavior of 
the system becomes more regular. 

Now we turn to the case when the atoms are prepared 
in the half-excited Dicke state \j, 0) at (for even values of 
N). The temporal behavior of the atomic inversion (J z ) 
is shown in Fig. |3| for n — 36 and N = 2,6, 10. Here 
we see a number of interesting phenomena. First of all, 
the amplitude of the Rabi oscillations is significantly sup- 
pressed. For n = 36, as N increases from 2 to 16, the 
relative amplitude Aq (the difference between the maxi- 
mum and minimum values of (J z ) in the very beginning 
of the evolution, divided by N) increases linearly from 
0.006 to 0.027, according to the empirical formula 



A = 0.0015(^ + 2). 



(3.6) 



For comparison, Aq is about 0.9 for to = j and about 
0.95 for to = —j, being almost independent of N. 
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FIG. 3. The atomic inversion (J z ) versus the scaled time 
r = gt for the initial coherent field state with n = 36 interact- 
ing with atoms prepared in the half-excited state \j, 0) a t: 
(a) AT = 2, (b) N = 6, (c) N = 10. 

Another important feature is the appearance of half- 
time revivals, i.e., for to = the revival time is actu- 
ally tr/2. In fact, tiny half-time revivals appear also for 
to = ±j, but their amplitudes are much smaller than the 
amplitudes of the first revivals at tr. As N increases 
from 2 to 16, the relative amplitude Ai/ 2 (the difference 
between the maximum and minimum values of (J z ) in 
the half-time revival, divided by N) increases monotoni- 
cally from 0.0018 to 0.0171 for m = j and from 0.0019 to 
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0.0354 for rn = —j. The ratio Ai/ 2 /Ai increases in the 
same range of N from 0.0033 to 0.0424 for m = j and 
from 0.0035 to 0.1677 for rn = —j. On the other hand, for 
rn = the half-time revival and the first revival are, for 
small values of N, of the same order of magnitude. For 
N = 2 we even find A1/2/A1 > 1. In contradistinction to 
the cases m = ±j, for m = the ratio A1/2/A1 decreases 
from 1.3427 to 0.0473 as N increases from 2 to 16. In this 
range of N the relative amplitude A 1 / 2 of the half-time 
revival increases slightly from 0.0037 to 0.0052, while the 
relative amplitude Ai of the first revival increases much 
faster: from 0.0028 to 0.1093. Starting from N = 4 the 
amplitude of the first revival exceeds that of the initial 
oscillations. The ratio Ai/Aq increases almost linearly 
from 0.4555 to 4.0476 as N increases from 2 to 16. The 
value (J z )c of the atomic inversion during the collapse 
is always positive for m = ±j. As N increases from 2 to 
16, {Jz)c increases from 0.0036 to 0.3569 for m = j and 
from 0.0035 to 0.5195 for m = —j. On the other hand, 
for m = the value of (J z )c is always negative and de- 
creases from —0.0070 to —0.2433 in the same range of N. 
We found that for m = the value of (J z )c can be we U 
approximated by the empirical formula 

( N \ z(Q) 

z(a) = &3 — exp[— (A^a + k$)], (3-7) 

where fa = 7.45, fc 2 = 11.16, fc 3 = 1-773, k 4 = 0.328, 
A15 = 1.681. We also see that in general the dynam- 
ics in the case m = is much more irregular than for 
rn = ±j. All these observations show that the structure 
of the phenomenon of collapses and revivals is essentially 
different for the half-excited Dicke state relative to the 
fully-excited state or the ground state. Inspecting the 
eigenvalues (see Table [jj), we see that the influence of 
the anharmonic corrections cannot explain this princi- 
pal difference. As expected, the behavior becomes more 
regular with increase of the mean photon number n. For 
odd numbers of atoms the half-excited Dicke state \j, 0) a t 
does not exist. However, for sufficiently large odd values 
of N the Dicke state \j, |) at exhibits properties very sim- 
ilar to those of the half-excited state. 

In order to explain the peculiar phenomena discussed 
above, we consider the exactly solvable case N = 2. For 
the fully-excited state (m = 1) we find 

( J z ) = J2 p n^r- \{n + 3) - (n + 1) cos(20 n ,ir) 

71 = n < X 

+8(n + l)(n + 2)cos(0„,ir)]. (3.8) 

Analogously, we obtain for the ground state (rn = — 1): 
00 ^ 

(J z ) = P ™ TP t( n ~ 2 ) ~ ncos(2fi n ,_iT) 

-8n(n - 1) cos(ft„ _it)] , (3.9) 



and for the half-excited Dicke state (rn = 0): 
00 ^ 

( J *> =-J2 P "7^2- I 1 - cos(2r!„, T)] . (3.10) 

Now it is clear why the oscillations are strongly sup- 
pressed for m = 0. Here both the constant and oscillatory 
terms have the prefactor P n j (4n+2), while for m — ±j in 
addition to the constant and oscillatory terms with pref- 
actors of the order P n /n there exists an oscillatory term 
with a prefactor of the order P n . For n = 36 the suppres- 
sion is by two orders of magnitude. For m = the oscil- 
latory term has the frequency 20 ni o = 2(n + 1/2) 1 / 2 , and 
the effective revival time is tr/2 = n(n+l/2) 1 / 2 . On the 
other hand, for m = ±j the leading oscillatory term with 
the prefactor of the order P n has the frequency f2 n ,±ij 
associated with the revival time tr, while the smaller os- 
cillatory term with the prefactor of the order P n /n has 
the double frequency, leading to the half-time revivals. 
This explains why for m — ±j the half-time revival is 
much smaller than the first revival while for m = 
both types of revivals are of the same order of magni- 
tude. The more irregular dynamics in the case m = 
can be explained by two reasons. First, the nonlinearity 
of the Rabi frequency (~ y/n) is less important for larger 
values of n, so the dynamics is more regular when the 
main contribution comes from larger n's. For m = ±j 
the prefactors of the leading terms are Poissonian, so for 
large enough n, the main contribution will come from 
the high-frequency terms, resulting in a regular behav- 
ior. However, for m = the prefactors are P n /(An + 2), 
so lower frequencies also contribute, which results in a 
less regular behavior. The second reason is that the re- 
vivals start to overlap much earlier if the revival time is 
halved. The same reason also leads to an additional ir- 
regularity in the case m — —j when the number of atoms 
is relatively large and the half-time revivals are not too 
small. 

The suppression of the revival amplitude for the 
half-excited Dicke state is similar to the trapping phe- 
nomenon which occurs for a single atom prepared in the 
equally weighted superposition state 2^ 1 / 2 (|+)±|— )) (see 
Rcf. |l4|]). In the latter case the population inversion is 
given by 

oc 

(J,) = -J2 P n[cOS 2 (VrT+lT) 

71 = 

+ ^. sin 2(V^TTr)]- \. (3.11) 
n + 1 2 

For large n, the Poissonian distribution is sharply peaked 
around n, and the two terms in the sum almost add up 
to 1. The remaining oscillating term has the prefactor of 
the order P n /n, so the amplitude of the Rabi oscillations 
is reduced by the factor of the order 1/n, As explained 
above, in such a situation the dynamics is less regular 
than for the case of initially unexcited or fully excited 
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atom. However, the revival time is not halted for the 
single atom in the equally weighted superposition state. 
In the single-atom case the population trapping occurs 
due to the destructive interference between the contribu- 
tions of the two levels, while for the two- atom half-excited 
Dicke state this phenomenon can be explained by the de- 
structive interference between the contributions of the 
two correlated atoms. 



IV. SQUEEZING OF THE RADIATION FIELD 

The coherent radiation field interacting with atoms 
can acquire interesting nonclassical properties such as 
sub-Poissonian photon statistics and squeezing. In the 
present paper we focus on the important quantum phe- 
nomenon of squeezing. The initially coherent cavity field 
can be squeezed when it interacts with a single atom 
|l5| , fL6|| , and squeezing in the revival-time regime can be 
very strong for large intensities of the field Jl7],[l8| . Butler 
and Drummond discovered [Q that in the Dicke model 
collective atomic effects can improve squeezing obtained 
for short interaction times, compared to the single-atom 
case. Squeezing in the Dicke model was also considered 
recently in a number of works pl[-p3[. Here we present a 
detailed study of squeezing of the initially coherent cavity 
field, comparing between the ground, fully-excited and 
half-excited initial atomic states. 

The field quadratures, q = (a^ + a)/V2 and p = 
i(a^ — a)/ \/2, satisfy the canonical commutation relation 
[q 7 p] = i- Then their uncertainties satisfy the Heisenberg 
relation: AqAp > 1/2, where (Aq) 2 = (q 2 ) — (q) 2 and 
similarly for (Ap) 2 . For the coherent field state (in par- 
ticular, for the vacuum state) the uncertainties are equal, 
Aqo = Apo = l/\/2, and an equality is achieved in the 
Heisenberg uncertainty relation. A field state is called 
squeezed, if the uncertainty of one of the quadratures is 
below the vacuum level, i.e., Aq < 1/V2 or Ap < 1/V2- 
Here we consider the squeezing parameter 

e = = V2Aq, (4.1) 

and search for £ < 1 that manifests field squeezing. 



A. Squeezing for the fully-excited and ground states 

1. The long-time behavior 

We first consider the temporal behavior of the squeez- 
ing parameter £ on the long-time scale, for the atoms 
prepared initially in the fully-excited and ground states 
(to = ±j). 




FIG. 4. The squeezing parameter £ versus the scaled time 
t — gt for the initial coherent field state with n = 36 in- 
teracting with N atoms (N = 2, 5, 10) prepared in (a) the 
fully-excited state |j, j) a t and (b) the ground state \j,—j)at- 

As demonstrated in Fig. |], this long-time behavior is 
fully correlated with the collapses and revivals of the Rabi 
oscillations. For even values of N, the squeezing param- 
eter £ oscillates achieving its minima at integer multiples 
of tr and maxima at half-odd multiples of tr. On the 
other hand, for odd values of N, we see that £ has minima 
at both half-odd and integer multiples of tr. This differ- 
ence between even and odd values of N can be explained 
by the fact that squeezing depends on two-photon transi- 
tions, when a pair of photons is simultaneously absorbed 
or emitted by a pair of atoms. The oscillations of £ decay 
with time and abandon their regular form. This decay 
is correlated with the loss of regularity in the behavior 
of the atomic inversion, caused by the overlaps of neigh- 
boring revivals. The behavior of £ is rather similar for 
to = j and to = —j, but in the latter case the decay of 
the oscillations and the irregularity occur earlier, as the 
revival time is shorter and the overlaps of revivals begin 
earlier. Also, for to = —j the decay rate of the squeezing 
oscillations increases rapidly with N, while for to = j the 
dependence of the decay rate on N is less pronounced. 
Inspecting the minima of £ in Fig. ^, we clearly see that 
their period becomes longer for m = j and shorter for 
to = — j a s N increases, in full accordance with the for- 
mula (3.5) for the revival time. 

For m = ±j, the minima of £ on the long-time scale are 
above the vacuum level, i.e., £ > 1. However, squeezing 
is achieved at short times, soon after the beginning of the 
interaction. In this region the squeezing parameter ex- 
hibits fast oscillations, with £ falling below 1. Therefore, 
for m — ±j we will focus on the behavior of squeezing 
on the short-time scale. 



G 



2. Squeezing on the short-time scale: The fully-excited state 

We first consider the case of the initial fully-excited 
atomic state (m = j). Squeezing is achieved for r < 2 
and appears at shorter times for larger values of N. 
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£ m (i.e., the maximum of squeezing), achieved during the 
time evolution, depends on N and n. 
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FIG. 7. The minimum value £ m of the squeezing parameter 
versus the coherent amplitude a for various values of N. 
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FIG. 5. The squeezing parameter f versus the scaled time 
t — gt for the initial coherent field state with n = 36 inter- 
acting with N atoms prepared in the fully-excited state: (a) 
N = 2, 5, 10, (b) N = 10, 60, 100. 

In Fig. | we see the short-time behavior of the squeez- 
ing parameter £ for n = 36 and various values of N. For 
relatively small values of N (N < 10), £ exhibits quite 
regular oscillations whose amplitude increases with N. 
However, for N ~ n and larger, the oscillations of £ be- 
come irregular. In Fig. ^| the short-time behavior of £ is 
shown for N = 14 and various values of n. The oscilla- 
tions of £ become more regular as n increases. 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
T 

FIG. 6. The squeezing parameter £ versus the scaled time 
t = gt for the initial coherent field state with n = 10, 15, 20 
interacting with 14 atoms prepared in the fully-excited state. 

It is interesting to investigate how the minimum value 




N 

FIG. 8. (a) The value « m jv of the coherent amplitude 
for which the squeezing parameter is minimized, versus iV 
numerical results (circles) and empirical fitting of Eq. (4.5) 



(line); (b) the corresponding value £ m jv of the squeezing pa- 
rameter versu s TV : numerical results (circles) and empirical 
fitting of Eq. (line). 

We first consider the dependence of £ m on a = \ffi for 
given N. As shown in Fig. ^, £ m has a minimum as a 
function of a, i.e., for given N there exists a value a m N 
for which the minimum value £ m j\r of the squeezing pa- 
rameter is achieved. Figure shows that a TO jv increases 
and { m jv decreases (squeezing improves) as N increases. 
We found that the dependence of a m M and £, m N on N 
at the considered range can be well approximated by the 
following empirical formulas: 



ot m N — a — be' 
£ m N=rN- s , 



cN 



(4.2) 
(4.3) 
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where a = 6.21, b = 4.03, c = 0.0471, r = 0.909, s 
0.156. 
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FIG. 9. The minimum value £ m of the squeezing parameter 
versus N for various values of a. 




FIG. 10. (a) The value N ma of the number of atoms, for 
which the squeezing parameter is minimized, versus a: nu- 
merical results (circles) and empirical fitting of Eq. (i.4) 
(line); (b) the corresponding value £ ma of the squeezing pa- 
rameter versus a: numerical results (circles) and empirical 
fitting of Eq. (O) (line). 



We next consider the dependence of £ m on N for given 
a. As shown in Fig. ^, £ m has a minimum as a func- 
tion of N, i.e., for given a there exists a value N ma for 
which the minimum value £ mQ of the squeezing parame- 
ter is achieved. Figure [lO] shows that N rna increases and 
£ma decreases (squeezing improves) as a increases. We 
found again that the dependence of N ma and £ mQ on a 
at the considered range can be well approximated by the 
following empirical formulas: 

N ma = -a + be ca , (4.4) 
£, ma =ra- s , (4.5) 

where a = 2.30, b = 1.216, c = 0.743, r = 1.159, s = 



0.490. Note that the relations ([h|) and ( p~4j ) are not the 
inverse of each other. This fact can be easily understood, 
if one imagine £ m as a two-dimensional function of N and 
a. When taking a section of £ m along the a axis (i.e., 
for a given N), one will find the minimum for a certain 
value of a. However, when fixing this value of a and 
going along the N axis, a minimum will be found, in 
general, for a different N. 

Finally, we would like to compare our numerical re- 
sults with approximate analytical expressions derived by 
Retamal et al. [p2| in the strong-field limit n~^> N. They 



found 1 22 1 the following expression for the minimum of 
the squeezing parameter, achieved during the time evo- 
lution, 



N 



N 



1/2 



(4.6) 



where a = l/\fe w 0.606 and 6 = 1 + a — a 2 — a 4 w 
1.103. Then, for given a, the squeezing parameter will 
be minimized by 



4a 





2.2a. 



(4.7) 
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FIG. 11. The minimum value £ m of the squeezing parame- 
ter versus iV for a — 6: (a) numerical result, (b) approximate 
analytical formula (4.6). 



In Fig. [IT] we compare our numerical results for £ m 
with the approximate formula ( |4.6| ) for a — 6 (n = 36). 
A good agreement is found only for very small values 
of N, while the values of N ma are absolutely different 
{N ma ~ 100 for our numerical calculations and N ma w 
13 for the analytical approximation). This discrepancy 
can be explained by the fact that the numerical values 
of N ma shown in Fig. [l(] do not satisfy the strong-field 
condition h ^> N. For the considered range of n (< 
36), the approximate solution ([0]) is not valid. From 
Eq. Q4.7D one can see that the strong-field condition will 
be satisfied for the optimal value N ma when N ma /n w 
2.2/cv <C 1. This means that Eq. ( [4.7[ ) gives a true value 
of N ma only for n ~ 10 4 or more. According to Eq. (fO|), 
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the absolute minimum of the squeezing parameter that 
can be achieved in the strong-field limit is £ « 0.58. For 
n = 36 (the maximum n we considered), the best value 
of squeezing is £ « 0.49 and it is achieved for N ma ss 100. 



3. Squeezing on the short-time scale: The ground state 

We next consider the case when the initial atomic state 
is the ground state (m = —])■ In this case we can dis- 
tinguish three different regimes: the strong-field regime 
(n S> N), the weak-field regime (n <C N), and the inter- 
mediate regime (n ~ N). 

The short-time behavior of the squeezing parameter 
£ in the strong-field regime is shown in Fig. [l2] (for 
a = 30 and N = 2,6, 16). We see that £ exhibits fast 
regular oscillations, whose frequency is almost indepen- 
dent of N (this is just the strong-field Rabi frequency 
fi^Lj = 2y/n- N/2 + 1/2, and the dependence on N is 
very weak because n is very large), but the amplitude in- 
creases with N. Therefore, the minimum value £ m of the 
squeezing parameter decreases as N increases. A similar 
behavior is found also for the case m = j in the strong- 
field regime. An interesting feature of this regime is that 
the minimum of £ is achieved after a relatively large num- 
ber of oscillations (the same time for different values of 
N) , while in the intermediate regime the minimum of £ is 
achieved, as a rule, in the first or second oscillation. As 
can be seen from Fig. [l3], in the strong-field regime the 
decrease of £ m is linear and can be well approximated by 
the empirical formula 



-f(a)N + k, 



(4.8) 



where the slope f(a) is a monotonically decreasing func- 
tion of a and the free term k w 1 is independent of a. 




N 

FIG. 13. The minimum value £ m of the squeezing parame- 
ter versus N for various values of a: numerical results (circles) 
and linear fitting of Eq. ( i.S ) (line). The fitting is good in the 
strong-field regime n> N. 








FIG. 12. The squeezing parameter £ versus the scaled time 
t — gt for the initial coherent field state with a = 30 inter- 
acting with N atoms (N — 2, 6, 16) prepared in the ground 
state. 



FIG. 14. The squeezing parameter £ versus the scaled time 
t — gt for the initial coherent field state with a = 1 interact- 
ing with N atoms prepared in the ground state: (a) JV = 10, 
(b) N = 20, (c) N = 30, (d) N = 40. 



It is well known J2jj2j,|3j,g3j that in the weak-field 
regime the behavior of the system is very regular (note 
that for the case of the fully-excited initial state the weak- 
field regime actually does not exist). In Fig. 14] we see the 
evolution of the squeezing parameter £ for relatively short 
times (for a — 1 and N — 10, 20, 30, 40). The oscillations 
of £ are quite regular and their frequency increases with 
N (this is just the weak-field Rabi frequency fi£"lj = 

2^J N — h/2 + 1/2), while the amplitude does not change 
significantly. The value of squeezing achieved in this limit 
is rather modest (as a rule, £ does not decrease below 
0.85). A perturbative analytical approach to the Dickc 
model dynamics in the weak-field regime was developed 
by Kozierowski and co-workers |2^j2^,^,^] and used for 
the study of squeezing in Ref. ]23fl . 
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FIG. 15. The minimum value £ m of the squeezing parame- 
ter versus N for various values of a. 

We focus our attention on properties of squeezing in 
the intermediate regime, where no analytical approxi- 
mation can be used. We study the minimum value £ m , 
achieved by the squeezing parameter during the time evo- 
lution, for N in the range between 2 and 40 and for a 
between 1.0 and 6.5. The results are presented in Fig. 
|l5| , where £ TO is plotted versus N for various values of a. 
The typical behavior for given a is as follows. Initially, 
£ m decreases slowly with N, but then steeply sinks down 
and acquires a minimum at a certain value of N. After 
the minimum, £ m slightly oscillates and then saturates 
for large iV's at an almost constant value. The region 
of linear decrease, occurring in the strong-field regime, 
appears only for a > 4. This leads to an additional mini- 
mum at small values of N, but it is less pronounced than 
the main minimum at larger iV's. 




FIG. 16. (a) The value N ma of the number of atoms, for 
which the squeezing parameter is minimized, versus a: nu- 
merical results (circles) and empirical fitting of Eq. (4.£) 
(line); (b) the corresponding value £ ma of the squeezing pa- 
rameter versu s a: numerical results (circles) and empirical 
fitting of Eq. (line). 



The value N ma , which gives optimal squeezing for 
given a, and the corresponding minimum value £ mct of 
the squeezing parameter are shown in Fig. |l6| as functions 
of a. As a increases, squeezing improves (£ mQ decreases) 
and the minimum occurs at larger N ma . The dependence 
of N ma and £ ma on a at the considered range can be well 
approximated by the following empirical formulas: 



Nina 



- ka 1 , 
x + ye 



(4.9) 
(4.10) 



where k = 1.382, 7 = 1.639, x = 0.476, y = 0.420, z = 
0.259. The largest value of a we considered is 6.5. If we 



assume that Eq. (4.10) remains valid for arbitrarily large 
a, then the absolute minimum of squeezing achievable 
with unexcited atoms is £ w 0.476. 



B. Squeezing for the half-excited state and the 
factorization approximation 

Finally, we discuss the case when the initial atomic 
state is the half-excited Dicke state (m — 0). The tem- 
poral behavior of the squeezing parameter £ on the long- 
time scale is shown in Fig. 0. We sec that there is a 
strong correlation between the collapses and revivals of 
the atomic inversion and the behavior of squeezing. 




FIG. 17. The squeezing parameter £ versus the scaled time 
t — gt for the initial coherent field state interacting with iV 
atoms prepared in the half-excited Dicke state: (a) n = 36, 
N = 2,4, 10, (b) n = 900, N = 2,4, 6, 10. 

For n = 36 and TV = 2, 4, we see in Fig. |l7|(a) that £ os- 
cillates with a period of order tr/2 (recall that for m = 
we find relatively strong half-time revivals). However, af- 
ter a number of periods, the amplitude of the oscillations 
decays and £ becomes nearly constant. This decay is cor- 
related with the loss of regularity in the behavior of the 
atomic inversion that happens when neighboring revivals 
overlap. For larger numbers of atoms, the behavior of the 
atomic inversion is very irregular, and so is the behavior 
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of squeezing. In particular, for N = 10, the oscillations 
of £ decay already after tr/2. The phenomenon of the 
collapses and revivals is quite regular for very large values 
of n. We see in Fig. 0(b) that for n = 900 the squeezing 
parameter also behaves very regularly. The structure of 
the oscillations of £ becomes more complicated (but keeps 
the regularity) as N increases. For all values of N, the 
squeezing parameter reaches deep minima at times just 
before integer multiples of tr/2. The larger the value of 
N, the closer is the minimum of £ to tr/2. (For N > 6, 
we also find additional minima of £, but they are not so 
deep.) 

The phenomenon of field squeezing for the half-excited 
Dickc state is drastically different from what we found for 
the fully-excited and ground states. In the former case 
(to = 0), for moderate values of n (a < 30), the min- 
ima of £ decrease below 1 (i.e., squeezing occurs) only 
for N = 2. This behavior is in contrast to the situa- 
tion in the latter case (to = where squeezing can 
be achieved for any N (with a proper choice of a) and 
is, moreover, enhanced by increasing N. However, just 
two atoms prepared in the half-excited Dicke state can 
produce quite strong squeezing. 

In Fig. |l| we see the temporal behavior of £ for N = 2 
and a = 6,10,16,30. The minima of £, which occur 
at times before integer multiples of tr/2, become deeper 
(i.e., squeezing improves) as n increases. It is very impor- 
tant to note that here squeezing is achieved at minima of 
the oscillations on the long-time scale. (Fast oscillations 
of £ on the short-time scale, which produce squeezing 
for m = ±j, are negligible for the case to = 0.) Conse- 
quently, for to = the duration of squeezing is essentially 
longer than for to = ±j. 



obtained for two atoms in the half-excited Dicke state 
has the same physical origin as revival-time squeezing 
in the Jaynes-Cummings model This long-time 

squeezing can be explained by means of the factorization 
approximation for the semiclassical atomic states p9| . 
This approximation is valid for a strong initial coherent 
field and for times short compared with To ~ n. 

In a semiclassical treatment, one replaces the boson op- 
erators of the field in the interaction Hamiltonian with 
c- numbers. The eigenstates \p) a t of this semiclassical 
Hamiltonian (called the semiclassical states) are just the 
eigenstates of the operator J + + J_ = 2 J x with the eigen- 
values X p = N — 2p (p = 0, 1, . . . , N). An atomic state 
(and, in particular, a Dicke state) can be expanded in 
the basis of the semiclassical states. The factorization 
approximation means that if the atoms are prepared ini- 
tially in a semiclassical state, then the total wave function 
of the system can be approximately written as a product 
of its field and atomic parts fl^ , |il"|| : 



l* P (*)> * \M*)U ® l*p(*))f- 



(4.11) 



Using the factorization approximation, one can ana- 
lytically estimate mean values of the field operators cor- 
responding to a specific initial semiclassical atomic state 
|p)at- In particular, one obtains (as usual, we assume 
that a = \ffi is real) : 



a exp 



iX r , 



2a 



1 4a 2 ) 



x exp < a 



i\ p T/4a i 



4a 3 



(4.12) 
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FIG. 18. The squeezing parameter £ versus the scaled time 
t — gt for the initial coherent field state with a = 6, 10, 16, 30 
interacting with two atoms prepared in the half-excited Dicke 
state: numerical results (solid line) and the factorization ap- 
proximation (dashed line). 

Actually, this kind of squeezing on the long-time scale 



/ 2\ 2 

{a )p = a exp 



x exp < a 



2a V 

i\ p T/2a 3 



iXpT 



2a 3 



(4.13) 



Also, (a^) p and (a t2 )p are given by complex conjugates 
of (a) p and (a 2 ) p , respectively, while {a^a) p = n can 
be taken constant. Then an approximate expression for 
the squeezing parameter can be easily obtained. For the 
semiclassical state |p) a t, one finds |2^] 



^ 2 = 1 



2n i 



^g-Tp/S _ e -T p 2 /16 

.3 f e -T p 2 /i6 _ 4e -T p 2 / 
16 V 



cos 



XpT 



+ 



2 

+2n ( 1 



-T 2 /16\ • / ^P T 

e J ? /iu sin 1 



-r 2 /i6 



(4.14) 



where T p = X p T/n. It can be seen that squeezing is 
achieved for any semiclassical state, except for those with 
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Ap = (then £ p = 1 is constant as long as the factoriza- 
tion approximation is valid). For A p ^ 0, the oscillations 
of the squeezing parameter achieve minima for times 

rW«0.9rW s 1.95r^ s ..., (4.15) 
where the revival time for the state |p) a t is given by 



_(z>) 



2-lTy/n T R 



I An 



I A. 



(4.16) 



The minima of £ p become deeper (i.e., squeezing im- 
proves) as n increases. 

Using formulas (4.12) and (4.13), one can calculate 



an approximate expression for the squeezing parameter 
for an initial atomic state which is a superposition of 
the semiclassical states. In particular, we will be inter- 
ested in approximate results for squeezing behavior of the 
Dicke states. First, note that in the Jaynes-Cummings 
model (N = 1) there are just two semiclassical states 
with A p = ±1, whose squeezing behavior is the same, 
for a good degree of accuracy. Consequently, any ini- 
tial atomic state in the Jaynes-Cummings model (and, in 
particular, the ground and excited states) will exhibit the 
same squeezing behavior (within the validity of the fac- 
torization approximation). These considerations explain 
the appearance of significant squeezing in the strong-field 
regime of the Jaynes-Cummings model for times near in- 
teger multiples of the revival time tr. 




a 

FIG. 19. The minimum value £ m of the squeezing parame- 
ter versus the coherent amplitude a for two atoms prepared 
in the half-excited Dicke state: numerical results (circles) and 
the factorization approximation (line). 

In the case of the Dicke model with two atoms, there 
are three semiclassical states; two of these states (with 
Ap = ±2) lead to squeezing, while for the third state 
(with Ap = 0) no squeezing is found (£ p = 1 is constant). 
The half-excited Dicke state |1,0) is a superposition of 
the two semiclassical states with A p = ±2, while the fully- 
excited and ground states also include the third semi- 



classical state with A p = 0. Consequently, the squeez- 
ing behavior of the half-excited Dicke state is very close 
to the behavior of a semiclassical state with |A P | = 2, 
while for the fully-excited and ground states squeezing 
is spoiled by the influence of the semiclassical state with 
Ap = 0. These considerations explain why squeezing on 
the long-time scale can be achieved only for the half- 
excited Dicke state 1 1 , 0) . The factorization approxima- 
tion also explains the appearance of the minima of the 
squeezing parameter at times close to integer multiples 
of tr/2, whose depth increases with n. 
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From Fig 

see that the factorization approximation describes very 
well the squeezing behavior of the two-atom half-excited 
Dicke state for large values of n. As shown in Fig. [H^, the 
minimum value £ m of the squeezing parameter decreases 
monotonically with a. For n > 10, the factorization ap- 
proximation is in excellent agreement with the numerical 
results. Similar results are also obtained in the strong- 
field regime of the Jaynes-Cummings model. For a = 38 
(the maximum value of a we considered) , the best value 
of the squeezing parameter is 0.3285. This is much better 
than the optimal values of squeezing which are obtained 
for the fully-excited or unexcited atoms when one takes 
N~10 2 . 

For atom numbers N > 3, any Dicke state will be 
a superposition of the semiclassical states with different 
values of |A P |. Since the minima of the squeezing parame- 
ter £p occur at different times for different |A P |, squeezing 
on the long-time scale will be spoiled for the Dicke states 
of three and more atoms. It can be shown f2^| , that for 
a Dicke state of TV-atom system the condition for the ex- 
istence of squeezing in the revival regime is n > (27V) 4 . 
Of course, the two-atom half-excited Dicke state is an ex- 
ception to this rule, as it involves only two semiclassical 
states with the same value of |A P |. However, the above 
condition is in good agreement with our numerical results 
for N > 3. 

As we see, the factorization approximation can be very 
useful for explaining many features of the field-atom in- 
teraction in the Dicke model. However, this approxima- 
tion fails to predict some interesting phenomena found 
by using numerical calculations. In particular, squeez- 
ing on the short-time scale, which clearly dominates for 
the fully-excited and ground atomic states (especially, for 
large values of N), is not predicted by the factorization 
approximation. We also observe that the suppression of 
the revival amplitude for the half-excited Dicke state can- 
not be described within this approximation. This can 
be re adily understood by recalling the exact expression 
( |3.10|) describing the evolution of the population inver- 
sion for the half-excited state |1,0). The amplitude of 
the Rabi oscillations here is of the order of 1/n, which is 
neglected in the factorization approximation. 
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V. CONCLUSIONS 

In this paper we considered in detail properties of the 
system of N two-level atoms interacting with a single- 
mode cavity field (the Dickc model). When the field is 
initially in the coherent state with a sufficiently large 
mean photon number, the dynamics of the system is quite 
regular, and the Rabi oscillations of the atomic inversion 
exhibit an interesting quantum phenomenon of collapses 
and revivals. We studied how this phenomenon is influ- 
enced by collective atomic effects. The main conclusion 
is that the role of the collective effects is determined by 
the initial atomic state. We found that by preparing 
just two atoms in the half-excited Dicke state one can 
cause greater effect on the behavior of the system than 
by collecting tens or even hundreds of excited or unex- 
cited atoms. In the phenomenon of the collapses and 
revivals, the half-excited Dicke state causes two basic ef- 
fects: the revival amplitude is strongly suppressed (anal- 
ogously to the trapping phenomenon for a single atom in 
the equally weighted superposition state) and the revival 
time is halved. The two-atom half-excited Dicke state 
also leads to very interesting squeezing behavior. It is the 
only Dicke state for N > 2 which exhibits strong squeez- 
ing on the long-time scale, similarly to the behavior found 
in the strong-field regime of the Jaynes-Cummings model. 
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TABLE I. The ratio / 5 i = hs/hi, that represents the 
anharmonicity of the Hamiltonian eigenvalues, for n = 36, 
N = 10, 20, 30 and m = ±j, 0. 
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5.0219 


20 


3 


5.0184 


30 




5.0159 


10 




5.0382 


20 


-3 
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5.0309 


30 




5.0360 
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